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Various properties of a material (viscosity, refractive index, etc.) 
can dramatically change after a phase transition. Investigating 
exactly when and how this change occurs potentially elucidates 
the dominant molecular forces working within the material, thus 
fueling significant scientific interest. Similarly, the analysis of 
idealized molecular models provides a framework for predicting 
the behavior of real, and often more complicated, materials. 
Understanding phase behavior is not just limited to fundamental 
research, as the ability to switch on or off various desirable 
properties with a simple "turn of the knob" can lead to interesting 
technological applications. 

Traditionally, intrinsic thermodynamic variables, such as 
temperature or pressure, act as the "knobs" that control phases. 
However, the application of external fields such as gravitational, 
magnetic, and electric fields, can play an equally vital roll since 
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they often offer a reversible, easily tunable, and localized method 
of control. The necessary key for practical use resides in a strong 
coupling between the field and a material property, for example 
the applied electric field with the material's dielectric response 
function. If such a strong coupling can be obtained, large changes 
to the phase diagram can occur even with a weak field. 

In this chapter, we will specifically focus on theoretical 
advancements for how electric fields induce a mixing-demixing 
phase transition between two dielectric fluids (oils), presenting 
results from both equilibrium and dynamics. Moreover, we 
will highlight how strong coupling between the field and 
dielectric function naturally occurs with nonuniform electric fields 
originating from small curved charged objects. 



1.1 Equilibrium phase behavior 



Using a mean-field approach, we consider a binary mixture of two 
fluids, A and B, in an electric field E, and write the total free 
energy T for a volume V as 



T = 



(fm + fe + fi)dV 



(1.1) 



where f m , f c , and fi are the free energy densities for mixing, 
electrostatics, and fluid-fluid interfaces, respectively. 

The fluids, in the absence of an electric field, can mix or 
demix due to a competition between entropy and enthalpy, where 
temperature T adjusts the relative balance. For concreteness, we 
consider the Bragg- Williams form 



Jm — 

VO 



</>ln 4> + (1 - (f,) ln(l -<f>) + xH 1 - 



(1.2) 



where k is Boltzmann's constant, vq is the molecular volume 
of both components, cf> such that < <f> < 1 is the volume 
fraction of component A , and x ~ i s the Flory interaction 
parameter Safran (1994)]. The first two terms account for entropy, 
while the third term accounts for the energy of mixing. 

Equation 11.21 gives rise to an upper critical solution 
temperature- type phase diagram. The shape of f m ((j>) transitions 
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from a single minimum at = 1/2 to a double minimum curve 
as T shifts from above to below the critical temperature T c . By 
using the well-known double tangent construction, we obtain the 
two binodal compositions 4*b for each temperature T < T c . These 
compositions mark the mixing-demixing boundary in the <f> — T 
plane. Specifically, fluids demix when is located between the two 
values of </>;,, or more simply stated "under the binodal curve". 
Finally, the phase diagram terminates at the mixture's critical 
point (<f, D , X c) = (1/2,2). 

In the following, we will employ a simplified form of eq. 11.21 
namely the Landau expansion of the mixing energy around <f) — <fic- 



kT 



1\ 2 4/ V 4 



(2-X) ( I +3 1 +const 



(1.3) 



Note that the quadratic term in the expansion changes sign at the 
critical value \c = 2. 

For electrostatics, the free energy density f c is given by 

/ o = ±i £Q £(0)|V^| 2 (1.4) 

where £o is the vacuum permittivity, s(<p) is the relative dielectric 
constant of the mixture, and ip is the electrostatic potential (E = 
—Vip). The positive (negative) sign corresponds to constant charge 
(potential) boundary conditions. The relation e{4>) can, in fact, 
be a complicated function. Using the simplest approximation, we 
assume a linear relation, e((f>) — {ea —£b)4> + £ b, where ea and Eb 
are the dielectric constants for pure fluids A and B, respectively. 

When considering the structure of the interface between phases, 
fi is required. This term accounts for the energ etic cost of 



composition gradients, and is given by Safranl (|1994I )]: 

/i = |^xA 2 |V0| 2 (1.5) 

where A is a constant characterizing the interface width. Notice 
that fi vanishes when the composition is uniform. 

To determine the equilibrium state in the presence of a field, we 
minimize T with respect to <f> and ip using calculus of variations 
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and obtain the following Euler-Lagrange equations 

S = V ■ [£ O e(0)W] - (1.6) 



ST _ kT 

8<j> v 

so de(» 



iv^r = M (i.7) 



2 ' 

The first equation is naturally Laplace's equation (Gauss's law) for 
the potential ip, while the second equation gives the composition 
distribution </>. Finally, e(4>) couples these two equations — in our 
case, de(0)/dr/> is constant. 

The Lagrange multiplier \x in eq. II .71 differentiates between open 
and closed systems. For a closed system (canonical ensemble), \i is 
adjusted to satisfy the mass conservation constraint: (r/>) = r/>o, 
where <po is the average composition. A closed system, whose 
volume increases to infinity, can be related to an open system 
in contact with a material reservoir (grand canonical ensemble). 
At this infinite size limit, the mass conservation constraint can 
be approximated as /i = fJ.o(4>o), basically the chemical potential 
that corresponds to the reservoir composition r/>o- The distinction 
between open and closed systems confers differences in phase 
behavior, as will be discussed below. 

Before continuing to nonuniform fields, we will briefly review 
changes in the phase diagram with uniform fields. The effect of a 
uniform electric field, Eo, on the mixture phase behav i or wa s first 
studied by Lan dau and Lifshi tz [Landau and Lifshitzj (|1957l )] and 
later by Onuki Onukil ([19951 )]. In the Landau theory, expansion 
of the electrostatic free energy leaves only a term proportional 
to {</> — <fi c ) 2 , which combines with the quadratic term in eq. 11.31 
to renormalize the critical temperature and the entire binodal 
curve. The theory predicts a critical temperature shift: AT C = 
vo£o£"Eg/(4fc) that is controlled by two free parameters Eo and 
e" = d 2 e((f) I 'd(j> 2 . If e" is greater (less) than zero, T c increases 
(decreases) and the electric field effect is that of demixing (mixing). 

Since vq is small in simple fluids, large electric fields are 
required to see an effect. For the typical maximal fields used in 
experiments (= 10 7 V/m), the predicted AT C is extraordinarily 
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small, on the order of miliKelvins. Experiments in low molecular 
weight binary mixtures agree with the theory on the magnitude of 
AT C , but yield conflicting r esults on the sign or direc tion of the 
shift |Debve and Klebothl (|l965l ): lOrzechowskil |l999l) ]. A more 
rigorous disc ussion about the effects of uniform electric fields is 
given in Ref. (Tsoril |2009h ]. 

Nonuniform electric fields, however, generate different results 
and can alter the mixing-demixing phase diagram c onsiderably. 



compa r ed to uniform fields of the same magnitude Tsori et al. 



( 2004 ): ! Marcus et al\ (jiool k ISamin and Tsoril |2009h ]. Large field 
gradients occur naturally in systems like microfluidic and nanoscale 
devices due to their small size and complex geometry. Detailed 
investigations of the phase transition have been conducted with 
three simple yet fundamental shapes — wedge, sphere, and cylinder. 
Analogous results occur between these shapes; therefore due to 
space constraints, we focus on a closed system consisting of two 
concentric cylinders with radii R\ and R2, where R2 — > 00 produces 
an open system. We impose cylindrical symmetry such that r/> = 
4>{r) and ip = ip(r), where r is the distance from the inner cylinder's 
center. Furthermore, the prescribed charge density a per unit area 
on the inner cylinder allows integration of Gauss's law to obtain 
an explicit expression for the electric field: E(r) — aRi/(eos(4>)r)r. 
Combining this result with E = —'Vtp in eq. 11.71 we obtain a single 
equation determining the composition profile </>(r): 



fi 



kT 



dfn 



de(»/d(/>. 
e(0) 2 1 



(1.8) 



where M = <r 2 v / (4kT c e ) is the dimensionless field, and f = r/R\ 
is the scaled distance. 

Inspection of eq. 11.81 shows two important differences from 
the case of uniform fields. First, the equation does not require a 
nonzero e" to shift the phase diagram. This holds irrespective of 
the three fundamental geometries. Recall that we, in fact, assumed 
a linear dependence for s((f>). Such a simple constitutive relation 
is insufficient for changing the phase diagram in a uniform electric 
field, resulting in AT C = 0. 

Second, the nonuniform electric field imposes a nonuniform 
"pull" on the fluid mixture, manifesting as an r-dependent total 
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Figure 1.1 The free energy density f(<j>, r) for a charged cylinder open 
system, (a) f(4>,r) — cfifj, versus <f> at distance r = n (dashed line), 
Ti (dash-doted line), and oo (solid line) for r/>o = 0.33, T/T c — 0.98, 
and M = 0.14. Symbols mark minima for each curve, (b) cj>(r) versus 
dimensionless distance f. Solid line is data from (a). Dash-dotted line 
has same r/>o and T as (a) but with M = 0.02. Dashed line has same 
0o and M as (a), but with T/T c — 0.995. (c) <j>{r) versus r for various 
values of M with <j> and T as in (a). 



free energy density f(4>,r) = f m + f c + f\. We will consider, for 
illustrative clarity, a fluid-fluid interface that is infinitely thin by 
specifically defining a vanishing interfacial term A = 0, which sets 
fi = 0. With this simplification, the behavior of / in an open 
system can be conceptualized as a competition between f m and f e . 
As r — > oo, the electric field is weak, f c — > 0, and f ~ f m governs 
fluid behavior. The solid line in Fig. 1.1(a) shows a typical example 
of f((f>,r) at a large value of r using </> = 0.33, T/T c — 0.98, and 
M = 0.14. The minimum of f(cj>, r), marked by a symbol, gives the 
value of 4>{r) as r —> oo, which in this case is 0.33. At the other 
distance extreme, r = Ri, the electric field is the strongest, and 
the dashed line in Fig. 1.1(a) shows the resulting f(<fi,r). Note the 
dramatic difference in the value of <f>(r) when the value of r is small 
(i?i) versus large. 

By finding the minima of / for all r, it is possible to construct 



the full concentration profile <fi(r), where the solid line in Fig. 1.1(b) 



corresponds to the data from Fig. 1.1(a) Whether or not a phase 
transition occurs resides in how the minimum of /(</>) changes as 
r varies between the two distance extremes. Specifically, if there 
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exists an r = r» where / contains two minima (see dash-dotted line 
in Fig. 1.1(a) I, then marks the interface between the two fluids. 



Figure 1.1(b) illustrates how the two minima in f(ri) translates 
into a discontinuity at <j){ri), thereby creating a distinct boundary 
between the two phases. 

Not all applied fields, however, induce phase separation. 
Figure 1.1(b) illustrates two such examples, where the dash-dotted 
line shows the same 4>q and T with a smaller M and the dashed 
line shows the same 4>q and M with a higher T. In both cases, 
/ contains a single minimum for all r, resulting in a smooth <p(r) 
profile. Despite the absence of a phase transition in these examples, 
the field still produces an effect, as the more polar (higher e) fluid 
accumulates near the high electric field. This phenomena can be 
thought of as the molecular version of the "dielectric rise" effect 
due to a dielectrophoretic force. 

Delving more deeply into the requirements for a phase 
transition, we vary M for a constant 4>o and T. Figure 1.1(c) shows 
that certain values of M induce a transition, whereas others do 
not. In fact, there exists a critical M c that marks the lowest M 
necessary for fluid-fluid separation. If an electric field can cause 
phase separation in a region of 4>o~T space above the binodal curve, 
a natural question arises: what is the new stability diagram for a 
particular M? This can be constructed by holding M constant and 
probing 4>a~T space for fluid-fluid demixing. Since the electric field 
breaks the symmetry of the free energy with respect to composition 
(4>o — > 1 — 4>o), the stability diagram is asymmetric with respect 
to (f>Q — <j) c . Figure 1.2(a) compares a typical stability curve for an 
open system, solid line, to the binodal curve, dashed line. Clearly, 
nonuniform fields can produce large changes to the phase diagram, 
in comparison with a uniform field. 

Figure 1.2(a) also includes the stability diagram for a typical 
closed system, dotted line, highlighting significant differences from 
open systems. Notably, the same M induces a weaker effect and 
produces a smaller stability diagram when 0o < 4> c , and the 
transition can occur when 0o > 4>c- These differences develop as 
a consequence of material conservation. Since there is no infinite 
bath from which to draw material, accumulation of <j) in high field 
regions leads to depletion of 4> m low field regions. Moreover, the 
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</>o M T/T c 

(a) (b) (c) 

Figure 1.2 Behavior of the fluid- fluid interface r^. (a) The stability 
curves in <f>o — T space for a constant M ~ 0.069 in an open (solid 
line) and closed (dotted line) cylindrical system. The binodal curve 
shown with dashes. (b,c) The location of n with respect to M with 
constant T/T c = 0.975 (b) and T with constant M = 0.02 (c) for 
open (open symbols) and closed (filled symbols) systems, where r/>o = 
0.29(v), 0.31(0), and 0.33(A). In (c), the dashed lines are the binodal 
temperatures for each value of r/>o, and the j/-axis is on a log scale. For 
all figures, R2/R1 = 5 in the closed system. 



penalty in / m grows faster than the energy gain in f c with changing 
4>o or T. Roughly speaking, however, the area under the stability 
curves for both open and closed systems increases (decreases) as 
M increases (decreases). 

Once an interface exists, several parameters control the location 
of n, for example </>o, T, M and i?2- In general, increases with 
increasing M (Fig. 1.1(c) and 1.2(b)), decreasing T (Fig. 1.2(c)), 
increasing R2 (not shown), and increasing r/>o (Fig. 1.2(b) I. The 



intere st ed reader can f ind m o re specific details i n Ref . [Tsori et al. 
|2004l k iMarcus et all (jiooj ): ISamin and Tsoril (|2009h ]. 



1.2 Phase separation dynamics 

Let's begin with a fluid-fluid mixture without an electric field and 
in equilibrium. When the electric field switches "on", the fluid 
mixture rearranges with time until it reaches a new equilibrium 
state. To describe the phase ordering dynamics, we follow the 
model H theoretical framework of the Hohenberg and Halperin 
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classification Hohenberg and Halperinl (|1977l )]. We assume that 
the electrostatic potential responds instantaneously to any changes 
in the composition, and supplement the model with Gauss's law to 
obtain 

V ■ (pv) = 



P 



- + (v.V)v 



V • [e e{<t>W\ = 



VP 



(1.9) 
(1.10) 
(1.11) 
(1.12) 



where v is the fluid hydrodynamic velocity field, P is the pressure, 
p is the fluid density, and r\ is fluid viscosity. The Cahn-Hilliard 
equation, eq. II. 9[ is the continuity equation for the mixture 
composition. The composition changes via two mechanisms: 1) 
a diffusive current that depends on gradients in the chemical 
potential /i = 8J-/S<p, where D is the diffusivity constant, and 
2) a convective current due to fluid velocity v. Equation 11.101 is 
the continuity equation for the fluid, while eq. II . 1 II is the Navier- 
Stokes equation t hat inc l udes a chemical potential-related body 
force -<jN5T/8<j) jOnukil (jjooj )]. Lastly, eq. EHis Gauss's law, 
which is again coupled through e(r/>) to the other equations. 

We will solve a simplified version of the model H dynamics, 
specifically in the over damped limit with no net fluid flow: 
v = 0. This is known as model B dynamics. The problem is 
now reduced to solving eq. 11.91 with only a diffusive current, 
and eq. 11.121 When T decreases from above to below the critical 
point in the absence of an electric field, diffusion governs the 
exchange of material, and the fluids phase separate locally into 
small domains that grow in time. Most notably, this process (at 
late times) can be described by a characteristic length such that 
the domain structure at all times is self-similar when rescaled by 
this length Cahn and Hilliard ( 19581 )]. Adding an electric field, as 
we shall see, alters this behavior. 

Returning to concentric cylinders in a closed system, we use 
the right-hand side of eq. 11.81 for the chemical potential in eq. 11.91 
and solve the radially symmetric problem with a dimensionless 
time t = {DkT / R\vo)t. When an electric field is turned on, 
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(c) 



Figure 1.3 Movement of r\ in time t. (a) (j>(r) versus dimensionless distance 
f at various dimensionless times t, where cf>0 = 0.33, T/T c = 0.965, M 0.31. 
(b) Symbols show interface location r, versus (for various M, where <f>0 = 0.33, 
T/T c = 0.965. Lines are fits to eg. 11.131 (c) The time constant t as a function of 
(j for various T/T c . For all data, R2/R1 = 5. 



material first accumulates near r « R±. If M > M c , a fluid- 
fluid interface emerges at r$(i,-) = Rt, travels outward to larger 
r, and asymptotically reaches the long-time steady-state location 



00) 



Figure 1.3(a) 



shows snapshots of 
mark n with time. 



(r) in time, 



while the symbols in Fig. 1.3(b) 

The behavior of can be approximated as an exponential 
relaxation 



i(*) = 



(n 



OexpI-Ct-tiVr] 



(1.13) 



which contains two free parameters, the time constant for 
relaxation r and the lag time for the interface to emerge U. The 



lines in Fig. 1.3(b) show fits to the data. In general, r decreases 



with increasing M (Fig. 1.3(c)), increasing T (Fig. 1.3(c)), and 
increasing <j) (not shown). 

In conclusion, we briefly reviewed important features of how 
nonuniform electric fields induce a fluid-fluid mixing-demixing 
phase transition. The advantage of nonuniform fields, over uniform 
fields, is that a phase transition can occur with only a simple 
dielectric difference between the fluids. Moreover, nonuniform fields 
can dramatically alter the transition temperature, compared with 
uniform fields. 

Because phase changes often lead to different material 
properties, the ability to easily control a phase transition often 
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translates into the ability to easily tune material behavior. For 
example, turning "on" an electric field can cause a homogeneous 
fluid-fluid mixture to phase separate, which creates a refractive 
index mismatch and produces an optical interface. Turning the field 
"off" reverses this process. Since nonuniform fields readily occur 
when electrical components are small, field-induced separation may 
have an important technological impact. 
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